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The complementarity of the liquid and plasma descriptions of the classical one-component plasma 
(OCP) is explored by studying wavevector and frequency dependent dynamical quantities: the 
dynamical structure factor (DSF), and the dynamic local field correction (LFC). Accurate Molecular 
Dynamics (MD) simulations are used to validate/test models of the DSF and LFC. Our simulations, 
which span the entire fluid regime (r = 0.1 — 175), show that the DSF is very well represented by a 
simple and well known memory function model of generalized hydrodynamics. On the other hand, 
the LFC, which we have computed using MD for the first time, is not well described by existing 
models. 
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I. INTRODUCTION 

The classical one-component plasma (OCP) is a stan- 
dard model in the study of strongly coupled plasmas, 
playing a conceptual role similar to that of the hard- 
sphere model in the theory of simple liquids. It is often 
used as a model of matter under extreme conditions, e.g. 
compact astrophysical objects. The OCP consists of a 
system of identical point charges Ze with mass m, inter- 
acting through the Coulomb potential, and immersed in 
a uniform background of opposite charge. In equilibrium, 
the system is characterised by the dimensionless coupling 
parameter T = (Ze) 2 /aksT, where a = (4-7rn/3) 1 ' 3 is 
the mean interparticle distance with n the particle den- 
sity and T the temperature. 

As r increases, the OCP changes from a nearly colli- 
sionless, gaseous regime for r < 1 through an increas- 
ingly correlated, dense fluid regime in which the system 
shares certain properties with ordinary liquids. In par- 
ticular, for T > 50 it has been found that the transport 
coefficients (diffusion, viscosity) of the OCP obey uni- 
versal laws satisfied by dense ordinary liquids Q. Other 
features of the OCP dynamics are not shared by ordi- 
nary liquids. Most notably, because of the long range 
Coulomb interactions, the system exhibits the character- 
istic behavior of plasmas: density imbalances lead to high 
frequency plasma oscillations, rather than low frequency 
sound waves. These high frequency plasma oscillations, 
not encountered in ordinary liquids, led Baus and Hansen 
to question the validity of the hydrodynamic limit of the 
OCP [J. In fact, it was recently shown that the hydrody- 
namic limit of the OCP is not applicable, even at large T 
values where high collisionality due to caging leads to liq- 
uidlike properties @- It is the fact that the OCP shares 
some, but not all, properties with ordinary liquids that 
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makes it a challenging yet fascinating system to study. 

In this paper we will explore the complementarity of 
the liquid and plasma descriptions of the OCP by study- 
ing the wavevector and frequency dependent dynamical 
structure factor (DSF), S(k,uj). The DSF contains com- 
plete information of the system dynamics at and near 
thermal equilibrium and is an important quantity be- 
cause of its connection to inelastic light and neutron scat- 
tering experiments [1, @. Two main approaches have 
been proposed for modeling the DSF in the fluid regime 
r < 175: the memory function approach and the dy- 
namic local field correction (LFC) approach. Largely due 
to the lack of 'exact' results (from numerical simulations) 
to compare to theoretical models of the memory function 
and LFC, it is not clear which of these approaches is more 
suitable for providing a description that is simple and ef- 
fective for a wide range of conditions. The purpose of 
this paper is to clarify this problem. 

The memory function approach - widely used for nor- 
mal liquids - represents a generalized hydrodynamics in 
which both equilibrium properties and transport coef- 
ficients that appear in the conventional hydrodynamic 
(Navier-Stokes) description are replaced by suitably de- 
fined wavevector and frequency dependent quantities. In 
this approach, the DSF is written in the form Q 

S(k,u) _ 1 <u>l > k 2 <P(k,oo) 

S{k) ~ ir [u 2 - <uj 2 > -uk 2 cf>"(k,uj)} 2 + [uk 2 <j>'(k,oj)] 2 ' 

(1) 

where S(k) is the static structure factor and < lo\ >= 
^ni~Wi?j' ^he quantities </>'(fc,cj) and c/)"(k,u)) are re- 
spectively the real and imaginary parts of the Laplace 
transform of the memory function <fi(k,t). In order to 
fully specify the DSF a model for the memory function 
is required. Of particular note is the Gaussian memory 
function model first applied to the OCP by Hansen et al. 
4] , which looked promising at the time of their study. 

The dynamics of the OCP can instead can be de- 
scribed in terms of the so-called dynamic local field cor- 
rection (LFC), G(k,L>). This approach is more common 
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to Coulomb systems e.g. the quantum electron gas [7|. 
The LFC is defined by its relation to the density response 
function of the system, x(k,uj) [8|, [9|, 



l-v(k)[l-G(k,uj)]xo(k,uj) 



(2) 



Here v(k) = A-K(Ze) 2 /k 2 is the Fourier transform of the 
Coulomb potential and \o(k,uj) is the density response 
function of an ideal gas, defined in Sec. IIII1 While the 
memory function is designed to extend the conventional 
hydrodynamic equations to finite wavevectors, the LFC 
is designed to correct the deficiencies of the mean field 
approximation (i.e. the Vlasov equation for the single 
particle distribution function, which describes the plasma 
oscillations but neglects any non-ideal or 'collisional' ef- 
fects). That is, setting G(k,ui) = gives the mean field 
approximation for the density response function; this 
gives a good description of the OCP dynamics in the 
weak coupling regime, r<l, only. A non-zero G(k,uj) 
represents correlation effects beyond the mean field ap- 
proximation. Models for the LFC have been proposed 
by Tanaka and Ichimaru [13] and by Hong and Kim [ll[ , 
but these models have barely been tested other than for 
a very few conditions in the original studies (and even 
for these conditions it was not clear how well the models 
agreed with the MD data). 

Since the density response function and DSF are re- 
lated through the fluctuation-dissipation theorem, 



ksT 

S(k,uj) = 3m{x(fc,w)} , 



(3) 



the LFC G(k, uj) is clearly related to the memory func- 
tion <j)(k,uj), albeit in a non-trivial way. In this paper 
we show that the memory function is a simpler quantity 
to model than the LFC. That is, a basic model for the 
memory function can describe both mean field and col- 
lisional effects that are characteristic of the DSF of the 
OCP, wheras an LFC that achieves this is much more 
complicated. Specifically, as shown in Sec. [Hi the Gaus- 
sian memory function model initially proposed by Hansen 
et al. reproduces the MD data for the DSF to remark- 
able accuracy across the entire fluid regime, and for all 
wavevectors k. In fact, the properties of the OCP mean 
that the model works even better than would be expected 
in the case of normal fluids. On the other hand, as shown 
in Sec. IIII1 the LFC has a more complex structure - 
for this reason it is not well described by the models 
mentioned previously. In order to reach these conclu- 
sions, we have performed highly accurate, large scale, 
state of the art molecular dynamics (MD) simulations 
for the intermediate scattering function F(k, t), and from 
this the dynamical structure factor, S(k,uj) [5], for a 
large number of T values spanning the entire fluid regime 
(0.1,0.3,1,5,8,9,9.5,10,11,50,120,160,175). We have used 
this new data to compute the LFC of the OCP with MD 
for the first time: calculation of G(k, uj) requires very ac- 
curate MD data which was not available before now. To 



conclude our study of OCP dynamics (Sec. IIV|) . we have 
extracted from our MD data the value of T at which 'neg- 
ative dispersion' of the OCP plasmon mode sets in; very 
recently there has been renewed interest in this particular 
aspect of OCP dynamics (T3 |. 



II. MEMORY FUNCTION MODEL 

The memory function expression in Eq. ([1} , which can 
be used to represent the DSF of any single component 
fluid, can be shown to be an exact result Q . In the case 
of the OCP, the ubiquitous plasmon peak in the DSF 
is ensured by the long wavelength (small k) behaviour 
of the < uj 2 > term in the denominator of Eq. (TT]); as 
k -)• 0, S(k) -> k 2 /k 2 D [l|, where k 2 D = 3T/a 2 , and hence 
< uj 2 >— > uip. This small k behaviour of S(k) is an 
essential distinction between OCP statics and those of 
an ordinary fluid - in the latter case, S(k) approaches 
the isothermal compressibility of the fluid in the limit 
k — > 0, which gives rise to a sound wave (rather than a 
plasma wave) at long wavelengths [l|. 

The memory function model first applied to the OCP 
by Hansen et al. Q consists of using the following Gaus- 
sian ansatz for the memory function, 

k 2 (f)(k,t) = fc 2 0(fc,O)exp(-7rt 2 /4T fc ) 

= [u>l{k)- < uj 2 >} exp(-^ 2 /4r fe ) , (4) 

where the initial value of the memory function is known 
exactly @ and uj^k) =< uj 4 > / < uj 2 > is given in 
terms of the frequency moments of S(k, uj) 



< uj'- >= 



uj n S(k, ui)dui . 



(5) 



Expressions for < uj° >, < uj 2 > and < uj 4 > in terms 
of the static structure factor S(k) and the radial distri- 
bution function g(r) for the OCP are given in the Ap- 
pendix. Here T k , appearing in Eq. is a wavevector 
dependent relaxation time. According to Eq. (Q}, the 
real and imaginary parts of the Laplace transform of the 
memory function are given by, respectively 0, fl3| , 



k 2 (f>'(k,uj) = [uJ 2 L (k)- < uj 2 >]T k e- T ^ 2 '* 



(6) 



and 



k 2 cj> [k,uj) 



2r k 



u 2 (k)- <uj 2 >}D(T k uj/V^), (7) 



where the Dawson function D{x) 
exp(-x 2 ) f*exp{y 2 )dy pjj]. 



A. Comparison between model and MD data 

The parameters < uj\ > and uj 2 L (k) that appear in the 
model can be obtained by computing S{k) (or equiva- 
lently g{r)) with MD and using the formulae given in the 
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Appendix for the frequency moments. The model then 
reduces to the determination of a single k dependent pa- 
rameter Tfc. The approach taken by Hansen et al. was to 
treat Tk as a parameter to be fitted to the MD spectrum 
oiS(k,u). 
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FIG. 1: (color online) MD results at selected T and ka val- 
ues (dots) and the result of the Gaussian model when the 
parameter Tk is fitted to the MD spectrum (dashed line). 



3 

to 

a. 
3 




0.035 
0.030 
0.025 
0.020 
0.015 
0.010 
0.005 
0.000 



_ r=i60 ; 


1 0.25 






fa»=1.85 j 


I 0.20 




fa=3.09 \ 


~ I 


1 015 








I 0.10 








t 0.05 












1 1 1 



0.0 



FIG. 2: (color online) MD results at selected F and ka values 
(dots) along with the result of the Gaussian model when only 
the parameter Tk is fitted to the MD spectrum (dashed line), 
and the result when < uf, >, < tdj j (k) > and Tk are all fitted 
(solid line). 

When we do this, we find in general that the model 
matches the MD data remarkably well for all T and k 
values (see Fig. [1]). In some cases, however, there are 
small discrepancies between the model and MD data, de- 
spite the fact that the model recreates the shape of the 
MD data very well (see Fig. [2j. Therefore, in order to 



determine whether these discrepancies are due to defi- 
ciencies in the model or inaccuracies in the parameters 
< ui 2 , > and < ijj\{k) > when computed with MD, we 
have separately fitted the model to the MD spectrum us- 
ing all three parameters < uj\ >, < u\(k) > and Tk. 
As shown in Fig. [21 this three parameter fit is an even 
better match to the MD data. Since the values of the 
parameters < lo 2 > and < w|(fc) > from the three pa- 
rameter fit agree very closely (within 10%) with those 
computed with MD, we conclude that the improvement 
in the agreement between the model when all three pa- 
rameters are fitted versus when only one is fitted is due 
to small inaccuracies when < uj 2 > and ui 2 L (k) are taken 
from the MD g(r) and S(k); the model is rather sensitive 
to the precise values of the frequency moments. That is, 
the one parameter fits are irrelevant as their comparison 
with the MD data is not indicative of the quality of the 
model. In Figs. [TU1 and [TTl we show only the model re- 
sults for when all three of the moments are used as fitting 
parameters. 

As shown in Figs. [TU] and [TTJ at sufficiently small fc, 
the MD data for S(k, w) exhibits a sharp plasmon peak at 
bj sa Lj p for all coupling strengths T. As k increases, the 
plasmon peak broadens until, at high k, S(k,u) reduces 
to a single central peak at uj = 0. The model accounts 
remarkably well for this entire evolution, particularly for 
T < 50 (Fig CGI. At higher values of T, the MD data 
does show some additional structure at intermediate k 
values (ka = 2.32 and 3.09) that the model cannot re- 
produce. For r > 120, a two peak structure is visible 
for ka = 2.32 and a three peak structure for ka — 3.09 
(Fig. ITTj) . The small high frequency peak for ka = 3.09 
is of particular interest - it does not appear to have been 
seen or commented upon in previous MD calculations. 
We believe that this peak is due to microscopic 'caging' 
effects. That is, at these lengthscales, the relatively high 
frequency oscillations of individual particles in the cages 
produced by their neighbors are imprinted on S(k,ui). 
This deduction is supported by previous work showing 
that for r > 100, a high frequency peak appears in the 
velocity autocorrelation function at ~ 0.9w p 0; this is 
exactly the position of the additional peak in S(k,ui). 
We note that although the model does not fully capture 
the additional structure in the MD data for these con- 
ditions, on average it does give a good account of the 
overall shape of the spectrum. 

As k increases further, SJk, ui) begins to reduce to its 
ideal gas limit S°(k,uj) 0, given by 



S°{k,u) 



{ m 
\2irk B Tk 2 



1/2 



exp 



2k B Tk' 



(8) 



As shown in Fig. [31 for T < 10, S(k,ui) is already close 
to S°(k,ui) for our highest k value (ka — 6.19). For 
r > 10, significant differences appear - these differences 
become greater as the coupling strength increases. This 
is to be expected, as at these higher coupling strengths 
oscillations the static structure factor S(k) persist well 
beyond ka = 6.19 (Fig. 0}. 
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FIG. 3: (color online) MD data for ka = 6.19 (dots) and the 
ideal gas limit given by Eq. ([H| (solid line). 




FIG. 4: (color online) The static structure factor S(k) for a 
range of F values. Clearly for T > 10 S(k) is still oscillating 
at our highest ka value. 



In any case, as shown in Figs. [TUlandlTTl the Gaussian 
model works well at our highest k value of ka = 6.19, 
regardless of whether or not this k value is sufficiently 
large for S(k, uj) to be close to its ideal gas limit. 



III. DYNAMIC LOCAL FIELD CORRECTION 

In the dynamic local field correction (LFC) approach, 
as given in Eq. @, the dynamics are described with 
reference to the density response function of an ideal gas, 



kuT 



2k B T k 



(9) 



where Z(x) = (1 — 2xD(x)) + iy/irx exp(— x 2 ) and D(x) 
is the Dawson function introduced in Sec. HU 

As mentioned previously, setting G(k,uj) = gives the 
mean field approximation for the density response func- 



tion; this only gives a good description of the OCP dy- 
namics in the weak coupling regime 1 < 1. A non-zero 
G(k,uj) represents correlation effects beyond the mean 
field approximation. One commonly used approximation 
is to replace the LFC by its uj = value G(k, 0) = G(k). 
The static local field correction G{k) is related to the 
static structure factor S(k) by 



G(k) = 1 + 



1 - 



1 



S(k) 



1 



ksT v(k) 



(10) 



An alternative scheme is to replace the dynamic local 
field correction by its high frequency limit 

G(k, oo) = G(k, uj -)• oo) = 21{k) , (11) 

where I{k), which depends on the radial distribution 
function g(r) of the OCP, is given in the Appendix. Re- 
placing G(k,uj) by either G(k) or G(fc,oo) results in a 
mean field approximation with an effective potential; it is 
well known that this type of scheme gives only a marginal 
improvement over the conventional mean field approxi- 
mation Thus, in order to describe well the dynamics 
of the OCP for r > 1, it is necessary to take collisions 
into account by having a frequency dependent dynamic 
local field correction. 

For the classical one-component plasma, two main for- 
mulations of the frequency dependent local field correc- 
tion have been given. The expression given by Tanaka 



and Ichimaru 10( , based on their viscoelastic formalism, 
interpolates between the known zero frequency and high 
frequency limits given in Eqs. (|10l) and (|11[) . 

G{k) — iuiTM{k)G{k, oo) 



G{k,w) 



(12) 



1 — iujTA[(k) 

In their prescription for computing the relaxation time 
T~M(k), Tanaka and Ichimaru considered either a Gaus- 
sian or Lorentzian ansatz [10(. In both of these cases, 
they used a kinetic equation to relate the shear viscos- 
ity to the local field correction; the unknown parameter 
appearing in rjvf(fc) was then chosen such that the es- 
timates of the shear viscosity available from MD at the 
time were matched as closely as possible (see 10] for fur- 
ther details). 

The other formulation, given by Hong and Kim 11], 
generates successive approximations for the LFC. The 
first order approximation is simply to replace G(k,uj) by 
G(k). The second order approximation is 

G{k, uj) = G(k) - \ \G(k) - G(k, oo)]Q (w^g^f 

(13) 

where Q(x) = l/Z(x) + 2x 2 - 1. Because Q(0) = 
and Q(x — > oo) = 2, like the model of Tanaka and Ichi- 
maru the Hong Kim model gives the correct zero and 
high frequency limits for the LFC. The third order ap- 
proximation involves the sixth moment of the dynamical 
structure factor, < oj e >; since this is difficult to compute 
theoretically, in the cases where the third order approx- 
imation was considered by Hong and Kim, < w 6 > was 
treated as an adjustable parameter [TT| . 
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A. Computing the Dynamic Local Field Correction 

In our MD simulations, we compute directly the inter- 
mediate scattering function F(k,t), The response func- 
tion is then given as 



x(M) = 



2k B T 



if t < 
dFML if t > o 



(14) 



dt 



which is simply the fluctuation dissipation theorem in the 
temporal domain (cf. Eq. ©)■ Numerically, we obtain 
the response function in the frequency domain, x(k,w), 
by taking the discrete Fourier transform of x(k,t). Fi- 
nally, we use the definition given in Eq. ^fy the compute 
the LFC. 

We find that the LFC is in general rather more difficult 
to compute with MD than the DSF; this is reflected in the 
less accurate and more noisy MD data we have obtained 
for G(k,ui). This is despite the fact that both S(k,u>) 
and G(k,u>) are derived from the same MD data for the 
intermediate scattering function F(k, t). In particular, it 
is difficult to obtain accurately the precise way in which 
the imaginary part of G(k, uS) decays to zero and the real 
part decays to its high frequency limit G(k, oo). This is 
because of the way in which G(k,uj) is defined (Eq. [2]): 
as to increases, both xo(k,u}) and x(k,uj) are small quan- 
tities. Despite these difficulties, our data is sufficiently 
good to allow for a comparison with the models of Tanaka 
and Ichimaru and Hong and Kim outlined above. 



B. MD results and comparison to models 

In Figs. [S] and El we have shown our MD results for 
the LFC at T = 10,50, 120 and 160 for small, intermedi- 
ate and large k respectively (ka = 1.02,3.09 and 6.19). 
In each of these figures we also show the model of Hong 
and Kim (Eq. (fTo) h Since the relaxation time Tjvf(fc) 
appearing in the model of Tanaka and Ichimaru (Eq. 
(fT2")0 was only given explicitly for V = 160 (see [JI|), 
we have shown their model - with both the Gaussian and 
the Lorentzian approximation for the relaxation time - 
for this coupling strength only 1 . 

At our small k value (ka = 1.02 - Fig. [5]), for all cou- 
pling strengths the model of Hong and Kim departs from 
G(k) and decays to its high frequency limit G(fc, oo) much 
faster than the MD data. At these long wavelengths 
(small fc), the OCP dynamics occur close to oj — uj p ; 
the local field correction describes the "collisional broad- 
ening" of the plasmon peak neglected in the mean field 
approximation (cf. Figs. [TU] and [TT]) . Therefore, the 
value of G(k,ui) should only be important for u> close 
to u) p . At our highest coupling strength of T = 160, 
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only) show G(k, oo). 



the model of Tanaka and Ichimaru works reasonably well 
for either a Gaussian or Lorentzian relaxation time. A 
reasonable estimate of the width of the plasmon peak is 
obtained, as noted previously [To| ■ 

At our intermediate k value (ka = 3.09 - Fig. [5]), 
G(k,w) shows rather more structure than at ka = 1.02, 
particularly at our largest coupling strengths of T = 120 
and 160. At these coupling strengths, the sharp vari- 
ation of both the real and imaginary parts of G(k,ui) 
around oj — uj p accounts for the high frequency peak in 
the dynamical structure factor S(k, uj) discussed in Sec 
III Al Again, for all T, the model of Hong and Kim de- 
parts from G(k) and decays to its high frequency limit 
G(k, oo) much faster than does the MD data. Further- 
more, at r = 160, the model of Tanaka and Ichimaru 
cannot capture the considerable structure in G(k,uj). 



1 In [l(J , the Gaussian approximation is referred to as "scheme I" , 
and the Lorentzian as "scheme II". 



At our large k value (ka = 6.19 - Fig. [7]), again none 
of the models seem to give a good description of G(k, ui). 
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FIG. 6: Same as Fig. [5] but for ka = 3.09. 



FIG. 7: Same as Fig. [5] but for ka = 6.19. 





r = 8 


r = 9 


r = 9.5 


r = 10 


r = 11 


b 


0.01948 


0.00887 


0.00221 


-0.00304 


-0.00523 


c 


-0.06467 


-0.06113 


-0.03124 


-0.02160 


-0.07313 


d 


-0.52751 


-0.44733 


-0.59515 


-0.61923 


-0.27276 





TABLE I: Parameters obtained from the fit for the plasmon 
peak position given in Eq. ()15|1 



IV. ONSET OF NEGATIVE DISPERSION 

Very recently, there has been renewed interest in the 
value of r at which 'negative dispersion' of the OCP 
plasmon mode sets in [P2J. Negative dispersion, refer- 
ing to duj(k)/dk < 0, where u) is the frequency and k the 
wavenumber of the plasmon mode, is a feature not pre- 
dicted by mean-field theory (i.e. the Vlasov equation for 
the single particle distribution function) . This anomalous 
effect, that represents an effect of the 'strong coupling' 
(r > 1) on the collective dynamics of the OCP, was first 
discovered in early computer simulations of the OCP [1] . 



3 




FIG. 8: (color online) Position of plasmon peak as obtained 
from MD simulations (open symbols) along with the least 
squares fits to the functional form in Eq. (|15[1 with the pa- 
rameters shown in Table U (solid lines) . 



A. MD results for plasmon peak position 

Figure [5] shows the plasmon peak position deter- 
mined from the MD simulations against ka for coupling 
strengths T = 8, 9, 9.5, 10, 11. As illustrated in Fig. [5J we 



obtained a good (least squares) fit of the MD results for 
w{k)/uj p to the polynomial 



w{k)/u p = 1 + ^{kaf + J(fca) 4 + ^{kaf . 



(15) 
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where higher order terms in ka were found to contribute 
negligibly to the quality of the fit. The values obtained 
for 6,c, and d at each coupling strength are given in Ta- 
ble HI From the b coefficient, we deduce that negative 
dispersion at long wavelengths sets in between T = 9.5 
and r = 10. Finally, in Fig. |9j we show the position of 
the plasmon peak for a larger range of T values. 
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FIG. 9: Position of plasmon peak as obtained from MD sim- 
ulations. 



and Ichimaru and Hong and Kim [ll[ - do not of- 
fer a good description for a wide range of conditions. As 
well as examining these two approaches, we used our MD 
data to accurately determine the coupling strength T at 
which the transition from positive to negative dispersion 
of the plasmon mode at long wavelengths takes place, as 
requested by a recently published study [l2[ . Aside from 
elucidating certain features of OCP dynamics, our MD 
results should find future application among practition- 
ers in the field of strongly coupled Coulomb systems. 



Appendix: Frequency moments of S(k, to) 



The wavevector dependent quantities, 



and 



<ujt >-- 



< or > 

<LU°> 



< 0T > 

< OJ 2 > 



(A.l) 



(A.2) 



are given in terms of the frequency moments of S(k,ui), 
defined as 



< w" >= 



uj n S(k, uj)duj . 



(A.3) 



V. CONCLUDING COMMENTS 

In this paper, we investigated two different approaches 
to describing the near equilibrium dynamical properties 
of the classical one-component plasma (OCP): the mem- 
ory function, which is a standard approach for normal 
liquids, and the dynamic local field correction (LFC), 
which is more familiar to plasma physics. Our study 
was centered around our highly accurate, state of the art 
Molecular Dynamics (MD) simulations for the interme- 
diate scattering function F(k,t). The accuracy of our 
MD data allowed us to compute not only the dynamical 
structure factor (DSF), S(k,uj), but also the dynamic lo- 
cal field correction (LFC), G(k,ui), the latter of which 
has to our knowledge never been computed before. We 
found that the memory function is rather more simple to 
model than the LFC: while the memory function is very 
well reproduced by a Gaussian for all coupling strengths 
r and wavevectors k, the LFC has considerably more 
structure. The more complex structure of the LFC is re- 
flected in the fact that current models - those of Tanaka 

2 Note that we define I(q) as in |fj|. This quantity has a numerical 
value of exactly half that defined in e.g. [8|], Eq. (3.37). 



The zeroth moment of S(k,u>) gives the static structure 
factor S(k) 



< uj° >= S(k) . 
The second moment is 



< ur > 



3r 



(A.4) 



(A.5) 



where q = ka with a = (3/(47rn)) 1 ^ 3 the average inter- 
particle spacing and uj p — y ^{Ze) 2 n/m is the plasma 
frequency. The fourth moment is [J] 



< oj 4 > 1 



sr 



q 2 - 2q 2 I(q) 



(A.6) 



with 



i(q)= / -=[<?(f)-i]x 



sin qr 3 cos qr 3 sin qr 



qr 



+ 



{qrf 



(qr)' 



dr . 



where r = r/a. 



(A.7) 
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